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Abstract 

We clarify the confusion in a recent paper by Aurenche, Zakharov and Zaraket (AZZ) 
[1] over the procedure and region of validity of the collinear approximation in the 
twist expansion approach to induced gluon emission in deeply inelastic scattering 
(DIS) off a nucleus target. We point out that, in this approach to the semi-inclusive 
spectrum, the transverse momentum £± of the induced gluon must be fixed in the 
collinear expansion in the transverse momentum k± of the initial partons, therefore 
the result is valid for (kj_) <C <C Q 2 - In the twist-four contribution, one can single 
out the double-hard term corresponding to collinear quark-gluon Compton scatter- 
ing which can be calculated independently of the collinear expansion approach. We 
will discuss the connection between the collinear approximation in the twist expan- 
sion approach and the small fey approximation of the results in the Light-Cone Path 
Integral (LCPI) and Gyulassy-Levai-Vitev (GLV) opacity expansion approach. We 
point out the misconstrued variable change by AZZ before the kr expansion in LCPI 
and opacity expansion approach, without which one obtains the same result for the 
induced gluon spectrum under collinear approximation as in the twist expansion 
approach. We also show that corrections beyond the collinear approximation to the 
transverse momentum integrated gluon spectrum within the static potential model 
in the GLV approach give rise to a logarithmic factor difference from the result of 
collinear approximation. 



1 Introduction 



Jet quenching due to parton energy loss in dense medium provides an excellent 
probe of the hot quark matter produced in high-energy heavy-ion collisions 
[2j. Such predicted phenomenon was indeed observed in experiments at the 
Relativistic Heavy- ion Collider (RHIC) [31H] . Phenomenological studies of the 
observed jet quenching depend crucially on our theoretical understanding of 
parton propagation and energy loss in dense medium which is dominated by 
induced gluon radiation via multiple scattering. There have been, therefore, a 
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Fig. 1. Diagrams of DIS with multiple gluon interaction 



plethora of theoretical studies of induced gluon emission from a propagating 
parton pioTITIIgll^irTU] . One approach to the study of induced gluon emission is 
based on higher-twist expansion of multiple parton scattering cross section in 
a nuclear medium in the framework of collinear factorization approximation 
[TUfTl] which can also be applied to parton propagation in hot medium. The 
advantage of higher-twist approach is the natural formulation of the parton 
propagation problem in terms of medium modified parton fragmentation func- 
tion which is the only physical observable of the jet quenching as a result of 
induced gluon emission and parton energy loss. 

The higher-twist approach to the problem of induced gluon emission is so far 
limited to twist-four contributions in the twist expansion. This is equivalent to 
the leading order approximation in the opacity expansion approach by Gyu- 
lassy, Levai and Vitev (GLV) [9]. The induced gluon spectra from these two 
approaches can be shown to be equivalent under the twist and opacity expan- 
sion approximations [12] ■ However, Aurenche, Zakharov and Zaraket (AZZ) 
claim in a recent note [1] that the higher-twist approach "fails" to produce 
the correct gluon spectra as in the Light-cone Path Integral (LCPI) and opac- 
ity expansion approach by Gyulassy, Levai and Vitev (GLV). We want to 
demonstrate that such an unjustified claim comes from their confusion over 
the procedure and validity of the collinear expansion in the higher-twist ap- 
proach. A misconstrued variable change in LCPI and GLV results also lead 
to their conclusion on the collinear approximation in the kr factorized formu- 
lation. We demonstrate that these results under the same approximation are 
equivalent to each other. 

For semi-inclusive cross section of induced gluon emission, the collinear expan- 
sion in terms of the initial parton transverse momentum kx must be made for 
a fixed value (kj) < <C Q 2 of the gluon's transverse momentum. For the 
total parton energy loss which involves integration over the gluon's transverse 
momentum, we will argue that the higher-twist expansion is a good approx- 
imation for large initial jet energy q~ ^> Ra^t), where Ra is the nuclear 
size. 
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2 Collinear factorization and twist expansion 



In the twist expansion approach to multiple parton scattering, one considers 
interaction between a fast parton and the target nucleus (or hot medium) via 
exchange of soft gluons as shown in Fig. [1] for the case of deeply inelastic 
scattering (DIS) off a large nucleus. In the framework of collinear factoriza- 
tion, one normally chooses a covariant gauge and considers A + as the largest 
component of the gauge field. Since the longitudinal momentum per target 
nucleon p = [p + , 0, 0j_] is the largest momentum component in the process, 
the twist expansion procedure involves expanding the hard part H n (kT) of 
the quark-n-gluon interaction in the transverse momentum kj< of the initial 
gluon fields and the transversely polarized gluon fields A±. Using generalized 
Ward identities one can relate the derivatives of the hard parts d^ T H{k T ) kT= Q 
of quark-longitudinal -gluon interaction and the collinear hard parts of quark- 
transverse-gluon interaction, and combine them to produce gauge invariant 
higher-twist contributions to the DIS cross section (see Ref. [13] for a detailed 
illustration). After integration over the transverse momentum of initial gluon 
fields, the final results are given by the convolution of collinear parton scatter- 
ing cross sections and transverse momentum integrated parton distribution or 
correlation matrix elements. The leading twist-four and nuclear enhanced con- 
tributions to the hadronic tensor of DIS off a nucleus involve two longitudinal 
gluon fields [H] as shown in Fig. [Tb and can be expressed as 



W2= I ^ d yi d y2^d 2 k T e-^Hr(k T ) 

x(A | ^0)J- A + (y^0 T )A + (y^^(y-) | A), (1) 

where we have suppressed the Lorentz indices of the electromagnetic currents 
and other kinematic variables in the perturbative hard part H 2 (fcr) of the 
multiple parton scattering with longitudinal gluon fields. Summations over 
color indices of the field operators in the matrix and average over the color 
indices of the initial state partons in the hard part are understood. In the 
collinear factorization scheme, one makes a collinear expansion of the hard 
part 

H-(k \ ?y-rm i p Hr{kT) i ,*44 ggTW ■ , (2] 
h 2 (k T )-H 2 (0 + ^-^—1^+— — ^-l^o + ---.(2) 

In order to clarify the confusion in Ref. [1] over the collinear expansion, it 
is important to emphasize here that the longitudinal gauge field A + is not a 
physical gluon field. Therefore, the hard part H 2 {k T ) does not correspond 
to quark interaction with physical gluons. This is apparent in the fact that 
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the collinear term in the above expansion does not vanish and is actually 
related to the hard part with no longitudinal gluon interaction (Fig. 
H 2 (0) = (—ig) 2 H . One can prove in general this is true for quark inter- 
action with any number of longitudinal gluons. After integration over the 
initial parton transverse momentum, which is another important part of the 
collinear factorization scheme, their contributions to the semi-inclusive cross 
section take the form 



H Q (A | ^(0) 7 - 



1 -ig I dy l A + (y l ) 




+ {-igf J dy^ J dy 2 A + {yl)A + {y 2 ) + ■ 




4>{y-) I A), 



(3) 



which becomes part of the leading twist contribution as the gauge link in the 
gauge invariant quark distribution function 

f%(x) = \J dy-e^ (A I ^(Oh+Pe-^cT ^ A+ ^(y-) | A) . (4) 



The interaction between a propagating quark and soft longitudinal gauge fields 
will only result in an eikonal line along the light-cone. This does not corre- 
spond to any physical scattering because one can get rid of it by choosing a 
proper (physical) gauge. This is the basic idea behind the proof of collinear 
factorization of the leading twist cross section of DIS and Drell-Yan processes 
by Collins, Soper and Sterman [I~6] . 

The contribution from the linear term in Eq. (j2J) of the collinear expansion 
vanishes for unpolarized targets. For the quadratic term of the expansion one 
can combine k\ with the longitudinal gauge fields A + A + and obtain a quark- 
gluon correlation distribution after partial integration over kx, 



J (2vr) 2 2tt yi y2 

x(A\mY k T A+ (y^°T)A + (y^^Hiy~)\A} 

f ^dy;dy;e^ + y- + ^ + ^^\A\m^F + *(y^ 

7T / dy^p A (y N )fl(x)x g G N (x g ), (5) 



where a factorized form of the quark-gluon correlation is assumed, Pa(vn) 
is the nucleon density distribution and Gn(x 9 ) the gluon distribution func- 
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tion per nucleon inside the nucleus. The momentum fraction x g carried by 
the initial gluon is determined by the kinematics of each individual process. 
Therefore, the leading twist-four contribution from the quadratic term of the 
collinear expansion to the semi-inclusive DIS cross section has a simple form 

7T J dy^p A (y N )^ d H Q 2} } kT ^ \k T =oJA(x)x g G N (x g ) (6) 



which has a simple and intuitive interpretation of partonic scattering between 
the fast quark and a physical collinear gluon since the contribution is propor- 
tional to the &/r-integrated gluon distribution inside the nucleus. The quadratic 
derivative term d\ T H 2 (k T ) kT=0 therefore corresponds to the hard part of the 
actual physical collinear quark-gluon scattering cross section. 

If one works in the physical gauge (A + = 0), the gauge link along the light- 
cone disappears (there will be transverse gauge link instead, see Ref. [17]). One 
therefore only has to consider quark interaction with the physical (transverse) 
gluons. In this case, the hard part iJ-f (&t), which will be convoluted with 
the (physical) gluon distribution, corresponds to partonic cross section of the 
physical quark-gluon scattering and is completely different from the hard part 
H 2 (kr) of quark interaction with longitudinal gluons in the covariant gauge. 
One can in fact prove in general the equivalence between the collinear hard 
part H^^O) in the physical gauge and the second derivative d\ T H 2 (£>r)fc T =o 
of the hard part in the covariant gauge [15] . We will illustrate this later in the 
case of induced gluon radiation. 



3 Induced gluon spectra in higher-twist expansion 



One can calculate nuclear modification to the dijet cross section in DIS [ISHS], 
differential direct photon [20] and Drell-Yan (DY) cross section [2~Tf22] in p+A 
collisions. The technique of higher-twist expansion has been also been applied 
to calculate the induced gluon spectrum due to multiple parton scattering 
in the DIS off a large nucleus or higher-twist contribution to gluon radiative 
correction to the semi- inclusive DIS cross section in Refs. pTj][TT] . In this case, 
one should keep the final gluon transverse momentum It fixed in the collinear 
expansion of the hard part in the initial parton transverse momentum kx- 
This is the origin of the flaw that leads to AZZ's conclusion in Ref. pQ about 
higher-twist approach to induced gluon spectrum. 

The dominant contribution in the twist expansion approach mainly comes 
from the process in Fig. [2j The hard part of the contribution to the semi- 
inclusive cross section before the collinear expansion is (Eq. (A. 18) in Ref. |llj). 
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I 

Fig. 2. Diagrams of induced gluon emission in DIS 



H — (kr) = 1 Qs r a 1 + — — c ™p + y-+i(*D+XL)p + (yi -y' 2 ) 



where 



XL 2p+g-z(l-*)' XD 2p+g-(l-*)' lj 



(4r - fc 



x D /z + x L = V " J _ '" Tj (9) 
2p + q z{± — z) 

Note that the overall phase factor quoted in Ref. [1] is different from exp[ixp + y~ 
i[xr> + xl)p + {%Ii — 1)2)] of the actual result. However, this does not affect the 
following argument. 



\2 



First, one notices that the collinear limit of the hard part H 2 (0) is finite 
and is the same as the vacuum gluon bremsstrahlung except the phase fac- 
tors. It combines with the collinear limits of hard parts from other cut dia- 
grams to form the gauge link in the quark distribution function in the vacuum 
gluon (leading twist) bremsstrahlung process in semi-inclusive DIS. To cal- 
culate higher-twist contributions in the twist expansion approach, one keeps 
the second order in the collinear expansion of the above hard parts in kx for 
fixed value of the gluon transverse momentum It- Therefore, one needs only to 
evaluate the second derivative of the above hard part, d\ T H 2 . The dominant 
contribution comes from differentiating the factor 1/(£t— &t) 2 - As explained in 
the previous section, one combines the quadratic term k\ and the longitudinal 
gluon fields to form the quark-gluon correlation function after partial 

integration. The corresponding contribution to the gluon radiation spectrum 
is 
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dN 2 a, i + (1-^)2 i 2™ 8 Tf g {x,x L ) 

dzdl% 2vr A z 1% N c fX(x) ' 1 J 

where 

T*(x,x L ) = J ^-dy^dy^il - e - iXLP+v *){l - e -^ P + {y--y^ e i{x+x L ) P +y- 
xe(-y 2 ")%--^^^|^(0) 7 + ^ + (l/ 2 )^( 2/ n^(2/1|A> ) (11) 

is the quark-gluon correlation function of the nucleus and f\{x) is the nuclear 
quark distribution function. 

The leading contribution to the second-order derivative d% H 2 from the 
phase factor in Eq. (J7J) will be linear in (yf — y2~)/q~ or y~ /q~ which in gen- 
eral are suppressed by a factor t\r N fq~ (r N is the nucleon size) or £^/xBP + q~ 
(xb is fractional momentum of the initial quark) relative to the above leading 
contribution for large jet energy q~. There are many other power-suppressed 
terms like these from other diagrams. They can be neglected for l\ <C Q 2 . 

Eq. (jTUI) has both hard-soft, double-hard scattering and their interferences. 
The double hard scattering is characterized by the finite momentum fraction 
xl carried by the initial gluon while in hard-soft contributions the initial 
gluon carries zero fractional momentum. This result is consistent with that in 
Ref [18] for nuclear enhancement of jet photoproduction where they consider 
only hard-soft and double hard scattering, but not their interferences which 
is not important for large values of l\. 

The higher-twist result also has a simple and intuitive partonic interpretation. 
One can assume the quark-gluon correlation has a factorized form (see Sec. II 
of Ref. [23] for details), 



Tq 9 {x, x L ) = A-k I dy N p A {vN) [f%{ x + x l) [xG n (x)] x=0 + f q N {x)x L G N (x L )} 
x [1 - cos(x L p + y^)}, (12) 

where f%(x) is the quark and G]y(x) the gluon distribution per nucleon, Pa(un) 
is the nucleon density at location y^ inside the nucleus A. 

The second term in the above factorized quark-gluon correlation corresponds 
to contribution from double hard scattering in which the initial gluon carries 
finite momentum fraction Xip + . The corresponding differential higher- twist 
gluon spectrum is then 
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dN. 



(H) 



dzd£\ 



dy N PA{yN 



ira 2 s C A 1 + (1 - z) 



£% n c 



da 



N 



d yNPA(VN)-^^[l - cos{x L p + y N )} 



x L G N (x L )[l - cos(x L p + y N )] 

(13) 



1 T 



which can be intuitively interpreted as the differential number of quark-gluon 
scattering in a collinear factorized form as the quark propagates inside the nu- 
cleus, where da^ g / ' dzd£\ is the collinear quark-gluon cross section on a nucleon 
target. 

One can derive the above double hard scattering contribution from the simple 
collinear factorized parton model by noting that Fig. [2] in this case is just 
the quark-gluon Compton scattering process. Considering a quark with mo- 
mentum q~ scattering with a gluon that carries a fractional momentum xp + , 
q(q) + g(xp) — > d(p') + g(£), the cross section can be written as 



o 4 d 3 £ 
da nh = ^\M\ 2 ab ^ cd (t/s, u/s)-——27rS[(xp + q - £) 2 } 



2s' 



(2vr)32£ 



\M\l b ^ cd (i/s,u/s)l-^-d£%5(l 



(4vr 



s 2 z{ 1 — z) 



X 



(14) 



where q = [0,q ,0] and xp = [xp + ,0,0] are momenta of the initial partons 
and 



02 
T 



2zq~ 



zq 



(15) 



is the momentum of the final gluon. With the given kinematics, the on-shell 
condition in the cross section can be recast as 



(xp + q — 



2(1 — z)xp + q ( 1 



XL 



x 



X L 



2z(l — z)p + q 



(16) 



The Mandelstam variables of the collision are, 



i 2 

s = (q + xp) 2 = 2xp + q~ = T — r, i = (£ — xp) 2 = —zs, 

z[l — z) 

u=(i-q) 2 = -(I - z )?± s = -(1 - z)s, (17) 
where the on-shell condition x = xl is used. 

With Eq. (fi~4l) and gluon distribution functions Gn(x), one can obtain the 
quark-gluon scattering contribution to the quark-nucleon cross section, 
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do-g g = J da gg G N (x)dx 



, XL G N {x L )\M\](i/s,u/s)-^^ ,d£ 2 



s 2 z(l — z) 



(18) 



Using the quark-gluon scattering matrix element 



|M| 



2 _ C A s 2 + u 2 C F s 2 + u 2 
qg ^g-Jf c —2 Jf c us 

C A l + {l-z) 2 CV (l + (l-z) 2 ) 
N c z 2 N c l-z 



(19) 



the quark-gluon cross section on a nucleon target is 



j N r> I \ 7Ta s 

da qg =X LGN{x L )-^r 



Ca m v . Of 2 

wS l - z) + w, z 



'l-z 



,2 



dzdt T , (20) 



z 

which is equivalent to the result in Eq. (|T3j) in the soft limit {z 0). 

The divergent factor in the induced gluon spectrum due to double hard 
scattering arises because of the collinear approximation in which we neglected 
the transverse momentum of the initial partons. This is related to the neglect 
of the /ct dependence of the phase factors in the hard parts in the collinear 
expansions and other higher-twist (larger than four) contributions. These ap- 
proximations are no longer valid at small values of i\ -C (k^) in the twist 
expansion approach, since many other contributions and processes will be- 
come important which are neglected in the above result. We have argued 
[TO] that contributions of these neglected terms could be approximated by 
substituting £^ — > l/£^(£^ + jj 2 ) and xlGn(xl) — ► (^l + x^)Gn{xl + x^) 
in the collinear result with x^ = fi 2 /2p + q~ and [i is the average transverse 
momentum of the medium gluon. However, the interference between double 
hard and soft rescattering processes suppresses the induced spectra for small 
£ 2 ± RA/2q~ ^ 1. Therefore, the final result in the collinear expansion will be 
a good approximation and insensitive to the regularization for large initial 
quark energy q~ ^> RA(k^). For corrections beyond the twist-four contribu- 
tion, one has to consider the nuclear broadening of the transverse momentum 
(k^) ~ RaQ {q is the jet transport parameter 



The double hard scattering process corresponds to elastic scattering in which 
there is a finite energy transfer (xl) from the medium gluon. With the leading 
order contribution to the medium gluon distribution xGn(x) ~ S(x — 1), the 
corresponding energy loss can be proved to be the same as the elastic energy 
loss [23]. This is a unique feature of the twist expansion approach that is not 
included in all other approaches (BDMPS and GLV). On the other hand, the 
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quantum non-locality of the quark-photon interaction in DIS has never been 
considered as an important feature of the twist expansion approach. 

One can similarly interpret the first term in Eq. (ITTj) which corresponds to the 
hard-soft process. Here the gluon radiation is induced by the hard photon- 
quark scattering and subsequently the radiated gluon interacts with a soft 
gluon from another nucleon with distribution xGn{x) x =o- One can compare 
this part of gluon spectrum 



dNP r _ ^a 2 s C A l + (l-z) 2 

dzW T = i dyNpAim) lfN- c z *G N (xU 

x[l - cos(x L p + y N )} (21) 

to the induced gluon spectra in LCPI and GLV approaches, especially after 
substitution xGn(x) x ~o — > x^G^x^) and £ T — > 1 / £^(£^ + fi 2 ) when the effect 
of the finite transverse momentum of the initial gluon is considered. 

Note that with higher order contributions to the medium gluon distribution 
function xGjy(x), the semi- inclusive spectrum from the double hard scatter- 
ing is similar to the hard-soft scattering with a correction on the order of 
(£t/Q 2 )[ x ^xGn(x)] x ~o which can be neglected for small values of £ T <C Q 2 . 



4 Induced gluon emission in A; T -factorized form 



The differential spectrum for induced gluon emission via interaction between 
the fast quark and the medium or initial gluon in the LCPI formulation is 
obtained in the factorized form, 



LCPI 



dzdi\ 



dy N pA{VN) 



d 2 k T 



xdG^{k T , x) 
din k\ 



H{kj 



'P T {k T -l T ) 



cos 



{l T - k T ) 2 y N 
2q-z{\ - z) 



(22) 



where q~ is the energy of the fast quark, fcy the transverse momentum of the 
medium or initial gluon, £t the transverse momentum of the emitted gluon (see 
Ref. [25|26] for similar formula in the GLV approach). Note that in Eqs. (12)- 
(14) in Ref. pQ the transverse momentum of the emitted gluon is replaced by 
£t — > £t ~ kr- Such replacement must be kept in mind when one expands the 
hard part in the initial transverse momentum k?- 



The above spectrum is proportional to the transverse momentum dependent 
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gluon distribution. One can interpret the corresponding hard part as the 
partonic cross section of quark scattering with a physical gluon. It vanishes 
H(kx) = at kx = 0, which is quite different from the hard part H 2 (kx) 
of quark and longitudinal gluon interaction in the higher-twist approach in 
a covariant gauge. Therefore, making a collinear expansion of this partonic 
cross section does NOT correspond to the collinear expansion of the hard part 
of quark and gluon interaction in the higher-twist approach in the covariant 
gauge. In order to compare to the results [Eq. ([6]) or (T2T1) ] in the collinear fac- 
torized formulation of higher-twist approach, one should integrate the above 
LCPI result over initial gluon's transverse momentum while keep the trans- 
verse momentum £t of the emitted gluon fixed. 

One can similarly make a small k? expansion of the hard part of the above 
LCPI result as in the higher-twist approach, 



H(k T ) = 4na 2 } kT ' ^ 



cos 



2q-z{\ 



+ 0(k T ). (23) 



After integrating over the initial gluon's transverse momentum and defining 
the /^-integrated gluon distribution function 



xG N (x) = I dk* xdGN J®' X \ (24) 



one can obtain exactly the same gluon spectrum as in Eq. (I2ip induced by 
hard-soft scattering in the higher-twist approach. The effect of higher order 
terms in the above expansion will be suppressed by powers of (k T )/£ T and 
(k T )£ T (Rji/ 'q~) 2 as in the collinear expansion in the higher-twist approach. 
Similarly, this approximation is no longer valid for small values of l\ < (k T ). 

To study the effect of these higher order terms in the kx expansion of the LCPI 
approach, one needs to know the form of /cr-dependent gluon distribution 
function Gx(k T , x). In the GLV approach [9], a static potential model was used 
for quark medium interaction. The induced gluon spectrum can be written as 



dN GLV _C A a s l + {l-zf 
dzd&r 



7T- 



dy N PA{yN. 



d 2 kr W 



vp ■ IT 



£ 2 T (k 7 



cos 



(4 + /i 2 ) 2 

{I T - k T ) 2 y~ N 
2q-z{\ - z) 



(25) 



where a q ^ is the quark-parton scattering cross section in the medium and fi 2 
is the screening mass in the static potential. 
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If one make a similar small kx expansion of the hard part in the above GLV 
result 



kx ■ It 



£ 2 T (k T -£ T ) 2 



cos 



(£ T - k T fy N 
2q-z{\ - z) 



cos 



-tUn 



2q-z{\ 



+ 0(k 



(26) 



the induced gluon spectra in Eq. (121!) from hard-soft scattering in the higher- 
twist approach can be recovered if one identifies the soft gluon distribution 
function as, 



2n 2 a t 



x^G N (x^ 



2l Q 2 

a qN /j log— 
/i 2 



(27) 



The above equation can be rewritten as, 



_ 1 X,G N { Xfl ) - J dk T 



x„G N (x„) = I dk 2 T -^fk 2 T . (28) 

T 



which relates the soft gluon distribution of the medium and the average trans- 
verse momentum weighted cross section in the higher- twist approach [BTTSB"] . 

Without the small kx approximation, one can make a variable change £' T = 
iT — kx in the integration over £' T , and complete the integration over the initial 
gluon transverse momentum kr, 



dN G LV C A a s 1 



dz 



2tx 



x 



zf 



dy N pA{yN)° q N^ 2 



d£' 2 



T 



1 — cos 



~tUn 



2q~z{\ - z) 



(29) 



The above result can be compared to the induced gluon spectrum in the higher- 
twist approach after regularization in £^, though one should note that £' T after 
the variable change in the above spectrum is no longer the true transverse 
momentum of the radiated gluon, which can be integrated over to obtain the 
transverse momentum integrated gluon spectrum dN/ dz. The above spectrum 
can be compared to the collinear factorized form of the spectrum in Eq. fl2Tj) . 
if one identifies 



2ttV 



X^Gn^X^) = (TqNlJ 2 



(30) 
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which differs from Eq. (I27p in the small kx approximation by a logarithmic 
factor 



Therefore, one can consider the above logarithmic factor as the conseqeunce 
of corrections beyond the small kx approximation within the static model of 
the GLV approach. 



5 Summary 

We have briefly reviewed in this paper the framework of twist expansion for 
the calculation of higher-twist contributions to semi-inclusive cross section of 
DIS off a nucleus target. We also illustrated the simple partonic picture of 
the higher-twist result for the induced gluon spectrum and demonstrated its 
equivalence to the quark-gluon scattering cross section in the collinear fac- 
torized framework. We pointed out AZZ's misunderstanding of the collinear 
expansion in this framework and their confusion over the procedure and region 
of validity of the collinear approximation. We showed that the small kx expan- 
sion in the LCPI approach without the misconstrued variable change leads to 
the same semi-inclusive gluon spectrum as in the higher-twist approach. We 
showed that corrections to the small kx (or collinear) approximation within 
the static potential model in the GLV approach lead to a logarithmic factor. 
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